Two Dimensional Ir-Cluster Lattices on Moire of Graphene with Ir(lll) 



Alpha N'Diaye/'Q Sebastian Bleikamp, 1 Peter J. Feibelman, 2 and Thomas Michely 1 

1 I. Physikalisches Institut, RWTH Aachen, 52056 Aachen, Germany 
2 Sandia National Laboratories, Albuquerque, NM 81185-1413, USA 
(Dated: February 6, 2008) 

Lattices of Ir clusters have been grown by vapor phase deposition on graphene moires on Ir(lll). 
The clusters are highly ordered, spatially and thermally stable below 500 K. Their narrow size 
distribution is tunable from 4 to about 130 atoms. A model for cluster binding to the graphene 
is presented based on scanning tunneling microscopy and density functional theory. The proposed 
binding mechanism suggests that similar cluster lattices might be grown of materials other than Ir. 

PACS numbers: 68.65.Cd, 81.16.Dn, 36.40. Sx, 61. 46. Be 



Fabrication of regular arrays of equally sized (monodis- 
perse) clusters on a flat substrate is a central goal for 
nanotechnology. Owing to their smallness, clusters differ 
from bulk materials in their chemical and physical prop- 
erties (cf. U). Because these properties depend strongly 
on size, monodisperse cluster arrays are optimal for fun- 
damental research and applications. Regular arrays of 
supported clusters are preferable to random ones, be- 
cause the identical environment of each cluster (e.g. dis- 
tances from their neighbors) produces a uniform response 
to external stimulation. Thus, in a regular array, one can 
use each cluster in the same way, either independently 
(e.g. for magnetic data storage), or by taking advantage 
of the coherent collective response of the array as a whole 
(e.g. in catalysis or for electrical transport). Recent ex- 
periments on two-dimensional regular arrays of clusters 
with a narrow size distribution have explored the size- 
dependent catalytic activity of Au-clusters |2|, the mag- 
netic properties of Co-clusters 3] and electrical transport 
through PbSe-clusters 4] , amply demonstrating the use- 
fulness of this approach. 

One route to cluster array fabrication is to deposit 
atoms or molecules from the vapor phase onto a "tem- 
plate," e.g., a substrate characterized by a periodic ar- 
ray of cluster nucleation sites, to which deposited par- 
ticles can diffuse. Examples are large unit cell super- 
structures of oxide films on metal single crystals p or 
regularly spaced steps and surface reconstructions Q. 
Here we demonstrate that graphene moires on an un- 
derlying dense-packed metal lattice act as templates for 
exceptionally well ordered cluster lattices with remark- 
able properties. 

Experiments were performed in an ultra high vacuum, 
variable temperature scanning tunneling microscopy 
(STM) apparatus with a base pressure in the 10 -11 mbar 
range. Sample cleaning was accomplished by cycles of 
flash annealing to 1500 K and sputtering by a mass sepa- 
rated 1.5 keV Xe + ion beam at 1100 K. Ethylene (5 L) ad- 
sorbed at room temperature was thermally decomposed 
at 1450 K resulting in the formation of large graphene 
flakes with sizes around 1000 A covering about 30% of 
the sample surface (compare [6]). By continuous expo- 




FIG. 1: (color online) (a) Atomic resolution STM topograph 
of graphene on Ir(lll). The rhombic moire unit cell is indi- 
cated by lines. Tunneling voltage applied to tip Ut = +0.2 V; 
tunneling current I t = 23 nA. (b) STM topograph after de- 
position of 0.02 ML Ir on graphene at 350 K; Ut = +0.2 V; 
It = 8nA (see text).(c) Schematic illustration of the DFT op- 
timized C(10xl0)/Ir(9x9) unit cell. Shading of the C atoms 
corresponds to their heights as calculated by DFT. 1st, 2nd, 
and 3rd layer Ir atoms are colored cyan, red and green. Hep- 
type region: full circle; fec-type region: short-dashed circle; 
atop- type region: dashed circle segments (see text). 



sure of the hot surface to ethylene, samples fully covered 
by graphene were also prepared. Ir was subsequently 
evaporated from a current heated Ir-wire, with a stan- 
dard deposition rate of 3.0 x 10 _3 ML/s, where 1 ML 
is the areal atomic density of the Ir(lll) surface. Pre- 
cise coverage calibration was performed by analysis of the 
fractional area of Ir islands in areas free of graphene. 

Fig. 1(a) is an STM topograph of the graphene moire 
on the underlying Ir(lll) substrate lattice. Pronounced 
bright regions are centered at the corners of its rhom- 



2 



bic unit cell, while the darker cell interior shows much 
less contrast variation. Superimposed on the moire 
is an array of small dark spots with the periodicity 
of the graphene lattice. Quantitative analysis of the 
moire by LEED and STM yields the following results: 
The dense-packed (lOlO)-directions of graphene and the 
unit cell vectors of the moire are parallel to the dense- 
packed (llO)-directions of Ir(lll) with angular scatters 
of (0 ± 0.26)° and (0 ± 2.6)°, respectively. The moire 
repeat vectors are of length, a m = (25.3 ± 0.5) A which 
equals 9.32 ±0.15 times the Ir nearest neighbor distance, 
ai r . The moire cell therefore accommodates A m = 87 ± 3 
Ir surface atoms. 

Fig. 1(b) displays Ir-clusters grown on graphene at 
350 K by Ir evaporation of 0.02 ML Ir. Note that the 
clusters are centered in down-pointing triangles of the 
moire's bright regions [compare Fig. 1(a)]. Some clus- 
ters were removed with the STM tip, prior to imaging, 
so we could view the corrugation of the moire together 
with clusters 7]. Owing to tip-surface interaction - of- 
ten during cluster removal - we also obtained images like 
Figs. 1(a) and 1(b) with inverted moire contrast. That 
is, pronounced dark regions surrounded the corners of 
the unit cell, while the brighter cell interior showed only 
weak brightness variation. In these images, not surpris- 
ingly, the clusters are found in down-pointing triangles 
of dark regions. 

To shed light on the C-Ir bonding, we optimized the 
geometry of a thin (3- or 4-layer) Ir(lll) slab with 
a graphene adlayer on its upper surface. The experi- 
mental graphene overlayer is not strictly commensurate 
with Ir(lll), but our model supercell, with a (10 x 10) 
graphene adlayer on a (9 x 9) Ir(lll) slab is an excellent 
approximation to reality. For the sake of interpreting 
the cluster bonding mechanism, we also conducted ex- 
ploratory calculations with 1, 3 and 4 Ir adatom clusters 
on the graphene/Ir(lll) supercell. 

We performed optimizations using the VASP, Den- 
sity Functional Theory (DFT) code jl Q in the PW91 
Generalized Gradient Approximation (GGA) with 
electron-core interactions represented by the Projector 
Augmented Wave approximation 0. The plane- 
wave cutoff was set to 400 eV, and Ir-Ir spacings in the 
rigid one or two lowest slab layers fixed at the PW91 
value for bulk Ir, 2.749 A. The surface brillouin zone was 
sampled by the point, f\ We accelerated electronic relax- 
ation by Fermi-level smearing (width = 0.2 eV) 0, and 
corrected for the contact potential difference associated 
with having a graphene adlayer on only one side of the Ir 
slab Systematic DFT error in lattice parameters is 
not a significant issue for the results. The PW91/GGA 
calculations imply that a graphene mesh need expand by 
< 0.4% to make a (10 x 10) graphene cell commensurate 
with a (9 x 9) Ir(lll) supercell. 

Graphene adsorbed in the experimental angular orien- 
tation presents three extremal regions for Ir cluster bond- 



ing, fee-, hep- and atop- type, named for whether an fec- 
or an hep-hollow, or an Ir atom shows through the local 
carbon hexagons. [Long-dashed, short-dashed and full 
circles circumscribe these regions in Fig. 1(c).] We 
infer that Ir clusters bind in hep-type regions based on 
the following logic: Fee- and hep- type regions are indis- 
tinguishable on a 1-layer Ir substrate, whereas atop re- 
gions, with all, rather than half the local C-atoms lying 
in Ir hollows, are structurally and electronically differ- 
ent. Adding a subsurface Ir layer is a second neighbor 
effect on the adsorbed C atoms, and thus perturbs their 
bonding weakly. Thus graphene layer properties in fec- 
and hep- type regions remain similar, and the atop- type 
region quite different. 

Accordingly, we identify the regions with the most pro- 
nounced contrast in the STM images [bright regions in 
Figs. 1(a) and 1(b)] as atop- type regions. Crystallogra- 
phy then implies that the clusters are adsorbed in the 
hep- type and not in the fee- type regions. Atomically re- 
solved topographs showing graphene and Ir(lll) side-by- 
side support this assignment: Assuming bright protru- 
sions on Ir-terraces to correspond to Ir atoms and dark 
spots on graphene to centers of carbon atom hexagons, 
by expanding the Ir lattice registry to the graphene one 
finds the bright regions in Figs. 1(a) and 1(b) again to 
be atop- type ones. It would add weight to our inferences 
to compare simulated STM images to experiment, but 
because of moire contrast sensitivity to tip condition, we 
have not. 

Lastly, worthy of mention is that for T < 160 K the fec- 
type regions of the unit cell also become partly populated 
with clusters indicating a second, shallower potential en- 
ergy minimum for Ir on graphene. Cluster adsorption at 
atop- type areas is never observed. 

Of considerable interest is the nature of the graphene 
bonds to the metal substrate. Given its strong sp 2 bond- 
ing, one expects that if graphene is to bind chemically 
to Ir (111), it must be as a result of C(2p 2 ) hybridization 
with the metal d-bands. Since the 2p 2 orbitals point 
along the normal to the graphene sheet, one expects C 
atoms in atop sites to form the strongest bonds to the 
metal through hybridization with Ir(d 32 2_ r 2) orbitals. 
The color coded C atom heights above the Ir(lll) sub- 
strate level in Fig. 1(c) bear out this expectation. The 
lowest-lying C atoms are in hep- and fee- type regions, 
3.77 A and 3.80 A above an underlying Ir, where many 
C atoms reside close to atop sites. The highest-lying C 
atom is found in the atop region 0.27 A higher than the 
lowest C of the graphene layer. The markedly different 
height of the C atoms in the atop-type region is in agree- 
ment with the experimental finding that the pronounced 
moire contrast extrema are atop-type regions. 

These results, incidentally, are found insensitive to 
whether the graphene layer sits atop a 3- or a 4-layer 
Ir(lll) slab, and whether, in the 4-layer case, only the 
bottom Ir layer is held rigid, or the bottom two. Adding 
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a fourth Ir layer has essentially no effect on the graphene 
layer's corrugation or binding energy, and changes its 
height above the Ir surface by no more than 0.05 A. 

The calculated average binding, 0.20 eV/C atom, rel- 
ative to a free graphene sheet and a clean metal slab, 
seems barely strong enough to indicate chemical bond 
formation. But this collective result likely masks forma- 
tion of relatively strong bonds in the hep- and fec-type 
regions of the moire cell, compensated by weak binding 
in atop- type areas. Evidence for this idea is the varia- 
tion, in the hep- type regions of Fig. 1, of the sizes of 
the openings through which red 2nd layer Ir atoms can 
be seen. That is a manifestation of in-plane shifts of 1st 
layer Ir atoms, to maximize hybridization with C(2p z ) 
orbit als. 

Structural optimization for Ir clusters has provided 
some insight into strong cluster bonding to hep-type and 
fee- type regions of the graphene moire. There, three out 
of six C-atoms in a carbon hexagon sit atop a substrate Ir 
and can form covalent bonds to it through hybridization 
of C 2p z orbitals with Ir(d 3 ^2_ r 2) orbitals (compare also 
[D|). This disturbs the graphene 7r-bonds, "activating" 
the remaining three C-atoms for bonding to the adclus- 
ter. Supporting this reasoning, cluster Ir's in the hep- 
and fec-type regions prefer binding not to C-atoms that 
lie atop Ir's, but to Ir's located over threefold substrate 
hollows. 

The sequence of STM topographs in Figs. 2(a)-2(e), 
taken after Ir deposition at 350 K together with the quan- 
titative analysis of the cluster density n and the aver- 
age cluster size in atoms s versus deposited amount 6 
in Fig. 2(f) suggest three regimes of cluster growth. In 
the nucleation regime represented by Fig. 2(a) s is only 
weakly dependent on with s ~ 4 — 5, while n in- 
creases nearly linearly with 0. Thus, at least monomers 
must be mobile on graphene at 350 K. The absence of 
monomers and dimers in the cluster size distribution for 
= 0.03 ML (not shown) indicates also dimer mobility 
at 350 K. The absence of a cluster denuded zone at the 
edges of the graphene flakes and the high n allow 
us to conclude that (i) adatoms and dimers are fairly 
well confined in the unit cell of their arrival or forma- 
tion and (ii) their intercell mobility is very low on the 
time scale of deposition (seconds). Therefore jumping of 
adatoms and dimers to neighboring cells continues after 
deposition during the « 15 minutes of slow cool down 
from 350 K to 300 K prior to STM imaging. After this 
time, as during imaging no more cluster mobility is ob- 
served, no adatoms and dimers are left. In the growth 
regime for 0.05 < < 1.50 ML framed by thin vertical 
dashed lines in Fig. 2(f) n is nearly [Fig. 2(b)] or exactly 
[Fig. 2(c)] equal to one cluster per moire cell and s oc 0. 
The growth regime is thus characterized by the slope one 
in the double logarithmic plot of s versus in Fig. 2(f). 
Compared to the Ir-Ir binding the Ir-C binding is weak. 
Thus, for energetic reasons at some size Ir clusters must 




FIG. 2: (color online) Cluster lattices on graphene flakes. 
The Ir-coverages are (a) 0.03 ML, (b) 0.10 ML, (c) 0.80 ML, 
(d) 1.50 ML and (e) 2.00 ML. Image height is always 550 A. 
Unfiltered data, (f) Cluster density n in % clusters per moire 
cell (full triangles) and average cluster size in atoms s versus 
deposited amount B in ML (full squares) . Note that s — A ™ @ . 
Lines to guide the eye. 



become three dimensional. Up to s = 6 clusters are pla- 
nar with an apparent height of about 2.2 A whereas for 
s « 25 two layer clusters are already as frequent as sin- 
gle layer ones. For = 0.80 ML [Fig. 2(c)] the average 
cluster height is 3.1 atomic layers and increases up 4.6 
atomic layers for = 1.50 ML [Fig. 2(d)]. The observa- 
tion of four to five layer high uncoalesced clusters with 
s = 130 and hexagonal top layer terrace suggests trun- 
cated polyhedra composed of the lowest surface energy 
{111}- and {100}-facets as typical cluster shape. Clus- 
ters in the growth regime are thermally stable up to 500 K 
with respect to intercell motion. Assuming an attempt 
frequency v$ ~ this translates to a potential energy 
minimum of a depth > 1.25 eV. Figure 2(d) with s = 130 
and = 1.50 ML also marks the transition to the coales- 
cence regime, already exhibiting a few coalesced clusters 
which extend over two graphene unit cells. It is apparent 
in Fig. 2(e) that coalesced clusters extending over sev- 
eral unit cells have a lower height, on average, compared 
to uncoalesced ones. The coalescence process causes a 
redistribution of material from upper to lower layers, re- 
sulting in a cluster height reduction [the average height 
in Fig. 2(e) is 3.7 layers]. Quantitatively the coalescence 
regime is identified in Fig. 2(f) by a decrease of n and 
by a superlinear increase of s. 
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FIG. 3: Experimental cluster size distribution (gray bars) and 
corresponding Poisson distribution (narrow black bars) after 
deposition of 0.80 ML at 350 K. Experimental size distribution 
adjusted to s as calculated from G. No other corrections for 
STM-tip effects. 



studies. It also might be possible to grow an oxide film 
on top of a cluster array, which then could be flaked off 
the graphene, resulting in an array of monodisperse tran- 
sition metal clusters in an oxide matrix. 
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Temperature variation alters cluster growth signifi- 
cantly. As mentioned for T < 160 K also the fcc-type 
regions are partly populated, while for T = 550 K and 
O = 0.80 ML five layer high, triangular clusters grow. 
These large clusters, with s = 650 atoms, extend over 
three hep-type regions. Although they are still in reg- 
istry with the moire, their positional order is worse than 
what is found in growth at 350 K. Extensive tempera- 
ture dependent measurements are underway to obtain 
quantitative information on the processes during cluster 
growth. 

A simple model for cluster formation would assume 
that all atoms deposited into a given unit cell form a 
single cluster growing in that cell, yielding a Poisson dis- 
tribution Pois(s) = fj-e -s of cluster sizes s. As discussed 
above, this model is inadequate in the nucleation regime 
for small s because of the intercell mobility of adatoms 
and dimers. However, since adatom and dimer intercell 
mobility is low during the comparatively short deposi- 
tion time, one expects the experimental distribution to 
approach Pois(s) as soon as O is large enough that the 
probability of having deposited fewer than three atoms 
into a cell is negligible. The experimental distribution is 
indeed fairly well described by Poig(«s) already for s = 9. 
We note that attachment of adatoms to clusters is irre- 
versible (consistent with the absence of cluster size 
changes during longtime STM imaging at 500 K). As for 
the Pois(s) the standard deviation is a = V^? the relative 
standard deviation cr r = | = 4| decreases monotonically 
with O. Fig. 3 exemplifies the agreement of the experi- 
mental cluster size distribution with the corresponding 
Poi7o(s) after deposition of 0.80 ML and demonstrates 
the narrow size distribution with cr r ,exp — 

12%. 

In conclusion, when graphene forms a moire on a 
transition metal substrate, with suitable unit cell size, 
one can locally functionalize the graphene with adclus- 
ters. Monodisperse cluster arrays on the relatively in- 
ert graphene surface open new opportunities for catalytic 
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